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Abstract 

Based on numerically accurate density functional theory (DFT) calculations, we systematically 
investigate the ground-state structure and the spin and orbital magnetism including the magnetic 
anisotropy energy (MAE) of 3d- and 4d-transition- metal dimer benzene complexes (TM2BZ, TM 
= Fe, Co, Ni, Ru, Rh, Pd; Bz = CgHg). These systems are chosen to model TM-dimer adsorption 
on graphene or on graphite. We find that Fe2, C02, N12, and Ru2 prefer the upright adsorption 
mode above the center of the benzene molecule, while RI12 and Pd2 are adsorbed parallel to the 
benzene plane. The ground state of C02BZ (with a dimer adsorption energy of about 1 eV) is 
well separated from other possible structures and spin states. In conjunction with similar results 
obtained by ab initio quantum chemical calculations, this implies that a stable C02BZ complex 
with Cq v symmetry is likely to exist. Chemical bonding to the carbon ring does not destroy the 
magnetic state and the characteristic level scheme of the cobalt dimer. Calculations including spin- 
orbit coupling show that the huge MAE of the free Co dimer is preserved in the C02BZ structure. 
The MAE predicted for this structure is much larger than the MAE of other magnetic molecules 
known hitherto, making it an interesting candidate for high-density magnetic recording. Among 
all the other investigated complexes, only RU2BZ shows a potential for strong-MAE applications, 
but it is not as stable as C02BZ. The electronic structure of the complexes is analyzed and the 
magnitude of their MAE is explained by perturbation theory. 

PACS numbers: 31.15.es, 75.30.Gw, 75.75.-c 
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I. INTRODUCTION 



Motivated by the ongoing quest for yet higher- density magnetic data storage in the con- 
text of the rapid advance of information technology, there is a continued search for nanoscopic 
magnetic structures with a large magnetic anisotropy energy (MAE) density by different 
experimental- and theoretical— methods. The MAE is the energy needed to turn the 
saturated magnetization of a system from one direction, usually the ground-state orientation, 
to another high-symmetry direction. Talking for simplicity about a nanoscopic system, we 
mean a canonical statistical ensemble of such systems, i.e. a macroscopically large number 
of identical non-interacting systems at equilibrium with a thermal bath. 

Thus, the MAE describes the stability of the magnetization direction against differently 
oriented external magnetic fields. In the simplest case of uniaxial anisotropy with an easy 
axis, the MAE is the energy barrier between two opposite equivalent directions of magne- 
tization with the lowest energy. A well-known example for this situation is bulk hep Co. 
Using the bi-stability of magnetic structures, a bit of information can be stored: for example, 
one stable direction of magnetization on a hard disc area may encode "0" , the other stable 
direction "1". It is generally accepted that long-term data storage requires that the total 
MAE of each magnetic particle should exceed 40 fcT,— where k is the Boltzmann constant 
and T is the temperature. 

Spin-orbit interaction in a magnetic state is the primary source of MAE.— The size of the 
spin-orbit coupling parameter of a given shell is a merely atomic property that depends on the 
atomic number. Hence, the spin-orbit splitting of atomic, molecular, or band states is fully 
determined by the character of these states in terms of atomic orbitals. On the other hand, 
the MAE as an energy difference depends sensitively on the particular electronic structure 
and, thus, on the geometry of the system. For instance, bulk hep Co shows a moderate 
MAE of 0.06 meV per atom. Larger MAE can be obtained in surface-supported structures. 
It can be considerably influenced by tuning the coordination and the hybridization through 
the choice of substrate and size of the deposited clusters. For instance, the deposition of 
single Co atoms on a Pt(lll) surface yields a record MAE of 9 meV per Co atom.- 

Further reduction of the dimensions leads into the realm of nano-particles and magnetic 
molecules. To give an example, Fe4 organometallic clusters with a propeller-like structure 
exhibit magnetic anisotropy barriers which can be tuned by altering the ligands and reach 
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up to 1.5 meV per cluster.- A large anisotropy barrier of 7 meV per cluster, generated by a 
deliberate structural distortion of the magnetic core with the help of bulky organic ligands, 
was recently reported for a complex within the Mn 6 family of single-molecule magnets.- 

Isolated magnetic dimers are the smallest chemical objects that possess a magnetic 
anisotropy as their energy depends on the relative orientation between dimer axis and mag- 
netic moment. Huge MAE values of up to 100 meV per atom were predicted for several 
transition-metal dimers (T12, Fe2, C02, N12, Zr2, TC2, RI12, Ir2, and Pt2).-^^^ However, 
it is impossible to utilize the huge MAE of dimers technologically unless they are bound to 
some medium. Our recent studies demonstrate that carbon-based substrates are suitable 
for this purpose.— Both benzene (Bz) and graphene are ideal support materials that do not 
spoil magnetism and the huge MAE of the Co dimer. The C02BZ complex has a ground-state 
structure with Cq v symmetry, in which the dimer is bound perpendicularly to the carbon 
plane. This hexagonal environment preserves the two-fold degenerate singly occupied high- 
est molecular orbital (HOMO) of the free Co dimer, which is responsible for the large dimer 
MAE.- As a result, C02BZ was predicted to show a magnetic anisotropy of the order of 100 
meV per molecule.— This finding may open a way to enhance the presently available area 
density of magnetic recording by 3 orders of magnitude. 

The present work has two aims. First, additional detailed results on C02BZ will be 



14. Second, related results for a whole 



presented that support the conclusions drawn in Ref. 
series of dimer-benzene complexes including the 3d and 4g? dimers Fe 2 , (Co 2 ,) Ni 2 , Ru 2 , Rh 2 , 
and Pd2 will be shown. Among these, only C02BZ and RU2BZ turn out to be interesting 
candidates for a potential strong-MAE application. 

Benzene-reacted metal dimers have been explored since the 1980s. Trevor et ai— studied 
the reaction products of benzene with gas-phase platinum clusters of different size. The 
existence of ES2BZ was verified with infrared spectroscopy by Ball et alM- after the reac- 
tion of iron atoms with cyclic hydrocarbons in an argon matrix. More recently, Bowen's 
group performed a series of mass spectrometry and photoelectron spectroscopic studies on 
iron-benzene,— iron- cor onene,— cobalt-benzene,— cobalt-pyrene,— cobalt-coronene,— and 
nickel-benzene^ cluster anions. The electron affinities and the vertical electron detachment 
energies were extracted from experimental spectra for the above complexes. By compar- 
ing measured spectra with results of density functional theory (DFT) calculations,— 1 ^ they 
proposed structure models for a part of these clusters. For example, a half-sandwich ground- 
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state structure {C§ v symmetry) was postulated for Fe 2 Bz~.— Their investigations suggest 
that carbon rings could be a suitable template to deposit small transition-metal clusters. 
The reaction of Rh+ cations with benzene was studied by Berg et alM These authors found 
that for n = 2 small amounts of Rh 2 Bz + were formed, while the main products were RhBz + 
+ Rh and Rh 2 C 6 H 4 + H 2 . Liittgens et al— measured the photoelectron detachment spectra 
of M 2 Bz~ (M = Pt, Pd, Pb) and resolved the electron affinities and ground-state vibra- 
tion energies of these complexes. By analyzing the vibration frequencies, they postulated 
a perpendicular arrangement of Pd and Pt dimers on Bz (Cq v symmetry) and a parallel 
coordination between Pb 2 and Bz (C 2 „ symmetry). 

On the theoretical side, DFT calculations were performed on several transition-metal ben- 
zene systems. Considering reaction products of iron atoms and benzene in low-temperature 
matrices, Parker recently calculated the energy of a number of isomers and simulated the 
related infrared spectra.— By comparison of calculated and measured spectra he concluded 
that Fe 2 Bz is formed at high iron concentrations. If the Fe dimer is assumed perpendicular 
to the benzene plane above the center of the carbon ring, the calculated infrared spectrum 
shows an excellent agreement with the experimental data, though the calculation finds a 
different ground-state isomer. 

The ground-state structures of Fe 2 -coronene,— Co 2 -pyrene,— and Co 2 -coronene2£ were 
also studied. For Fe 2 -coronene, a ground state with a total spin 5 = 3 was found with 
three quasi-degenerate isomers, where the Fe dimer is oriented either parallel or perpendic- 
ular to the coronene plane.— For both Co 2 -pyrene and Co 2 -coronene, the ground state was 
found to be 5 = 2 with perpendicular orientation of the Co dimer at a position above a 
peripheral C-C bridge.— 1 ^ Earlier work in this field is due to Senapati et al.— who studied 
neutral and cationic Fe 2 -coronene complexes but discussed only parallel adsorption modes. 
Also, Rao and Jena studied the geometries and magnetic moments of neutral and ionic 
Ni n Bz m complexes.— 1 ^ They predicted a parallel adsorption mode between the Ni dimer 
and benzene, but it is unclear whether a perpendicular geometry was considered or not. 
The interaction of benzene with Rh + and with Rh^ was investigated by Majumdar et al.— 
For the physisorbed dimer cation, Rh 2 Bz + , the minimum energy geometry has C 2 „ symme- 
try with the two rhodium atoms lying horizontally above the benzene at the C-C bridge 
sites. 

Further research activities were devoted to transition-metal dimers interacting with 
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graphene or with fullerenes. Interaction of silver and gold adatoms and dimers with graphite 
or graphene was studied,— 1 ^ and a perpendicular orientation of gold dimers on graphene 
was predicted.— Also, both structure and spin magnetic properties of 3c? transition-metal 
adatoms and dimers on graphite were investigated by Duffy et al.,— but a possible per- 
pendicular arrangement of the dimers was not considered. DFT calculations for palladium 
clusters supported on graphene^ and on find that the two atoms of an adsorbed Pd 

dimer are located on bridge sites, i.e., on top of C-C bonds. Recently, two detailed theoret- 
ical studies of Fe, Co, and Ni adatoms and dimers adsorbed on graphene were published: 
Johll et al— found that the most stable structure for all considered dimers, Fe2, C02, and 
Ni 2 , has a dimer axis oriented perpendicularly to the graphene plane and placed at the hole 
site. An enhancement of the magnetic moment for the atom farther from the graphene 
was predicted,— compared with the free dimer. Cao et a/.— found the same ground-state 
geometry, if the density-gradient corrected functional according to Perdew, Burke, and Ernz- 
erhof— was used, but they note that partly different results were obtained by using the local 
spin-density approximation (LSDA). 

DFT calculations do not only allow to predict the ground-state geometry and spin of a 
magnetic complex, but may also provide a basic understanding of the electronic structure 
and of the orbital magnetic properties like orbital magnetic moment and MAE. We are 
however not aware of any published calculations of the magnetic anisotropy of transition- 



metal dimers on carbon-based systems except our recent letter, Ref. |14|. In the following, 
we will demonstrate by DFT calculations that carbon hexagons are suitable hosts, where 
adsorbed transition-metal dimers may preserve their exceptional magnetic anisotropy. 

The investigated TM 2 Bz complexes are meant to serve as model structures for the ad- 
sorption of transition-metal dimers on the surface of graphite, on graphene, or on other 
carbon structures including molecular systems. For this reason, we only consider so-called 
physisorption (adsorption without expelling other atoms, e.g., hydrogen) as opposed to 
chemisorption that includes the possibility of de-hydrogenation.— Most probably, the pre- 
sented predictions can only be verified under ultra-high vacuum conditions. Any interaction 
of the transition- metal atoms with, e.g., oxygen may deteriorate the specific structure and 
the related magnetic state we are focusing on. 

The paper is organized as follows. In Sec. [Til the calculation method and computational 
details are explained. Sec. [TTT] compiles all results and related discussion: structure opti- 
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mization and stability of the ground states, analysis in terms of the bonding mechanism, 
the spin and orbital moments, and the strength of the MAE. The origin of the huge MAE 
in some of these molecules is also explained. Finally, the paper is summarized in Sec. IV. 
The appendix contains a description of auxiliary calculations. 

II. METHOD AND COMPUTATIONAL DETAILS 

The DFT calculations were performed with an highly accurate all-electron full-potential 
local-orbital scheme (FPLO),— release 8.00-31.— The code is based on a linear combination 
of overlapping nonorthogonal orbitals with a compact support. The molecular mode of 
FPLO with free boundary conditions was used. The presented data were obtained using 
the generalized gradient approximation (GGA) with a parameterized exchange-correlation 
functional according to Perdew, Burke, and Ernzerhof.— All results were checked against 
additional calculations using the LSDA in the parameterization by Perdew and Wang.— 
The dimer adsorption energy calculated by GGA is in all cases about 1 eV smaller than 
the related LSDA energy, but both approaches find the same ground-state structure type. 
Also, the ground-state spin obtained with GGA or LSDA is the same for all systems except 
Rh 2 Bz, where LSDA yields a non-magnetic ground state and GGA yields an S — 1 ground 
state. In both cases, however, the energies of the 5 = and 5 = 1 states are very close. 

The molecular levels were occupied according to a Fermi-Dirac distribution in order to 
ensure the convergence of the Kohn-Sham equations. The presented results were obtained 
with a broadening temperature of T = 100 K. The basis set comprised 3<i-transition-metal 
(3s, 3p, 3d, 4s, 4p, 4<i, 5s), 4o?-transit ion- metal (4s, 4p, 4d, 5s, 5p, 5d, 6s), carbon (Is, 2s, 
2p, 3s, 3p, 3d), and hydrogen (Is, 2s, 2p) states. Lower-lying states of the transition-metal 
atoms were treated as core states. 

Geometry optimization was carried out with a scalar relativistic scheme. To find out 
the lowest-energy geometry and spin magnetic state of each TM 2 Bz complex, three possible 
high-symmetry structures (Fig. [1]) were optimized for S = 0, 1, 2, and 3 (total spin moment, 
fis = 0, 2, 4, and 6/ie), and for different initial spin arrangements (ferro- and ferrimagnetic) . 
The point group symmetry Cq v was applied for the configuration shown in Fig.[TJ(a), while for 
the structures depicted in Figs. H(b) and[TJ(c), Ci v wa s used for the structure optimization. 

Previous theoretical investigations of metal-benzene systems have shown that the struc- 
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tural changes of the benzene plane due to the metal-benzene interaction are negligible.- 1 ^— 
Thus, in nearly all our calculations the positions of the C and H atoms were fixed with C-C 
bond length 1.40 A and C-H bond length 1.09 A. Exceptions from this strategy are reported 
below. 

The C02BZ complex is of particular interest.— To make sure that the correct ground state 
was found for this system, we utilized the pseudopotential code ESPRESSO-4.0.1— and 
cross-checked the above calculations by full optimization of all the atomic positions starting 
from 14 kinds of initial structures. The computational details and a brief description of the 
results are given in the appendix. 

A quantity used to judge the stability of the considered structures is the adsorption energy 
(E a d) for a dimer entity attached to a benzene molecule, which is defined as 

E ad = E tot {Bz) + E tot (TM 2 ) - E tot (TM 2 Bz) . (1) 

Here E to t refers to the respective total energy of the species indicated in the parentheses. 
Negative values of E a d mean that the TM2BZ complex is unstable. 

To evaluate the orbital magnetic moment and the MAE, spin-orbit coupling has to be 
included in the calculation. However, standard (quasi-)local DFT approximations like LSDA 
or GGA do not include orbital-dependent exchange effects.— Thus, the orbital moments and 
the MAE are usually underestimated by these approaches. The orbital polarization (OP) 
correction^ is a frequently applied method to cure this problem. As a matter of experience, 
the MAE evaluated with standard LSDA or GGA approximation gives a lower estimate to 
the expected MAE, while the value obtained by including the OP correction provides an 



upper estimate. Experimental values of the MAE are most probably located 
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Detween these 
and 3 and also for 



lower and upper estimates. This has been demonstrated, e.g., in Refs 
the special case of Co atoms in different chemical and structural surroundings in Ref. [14, 
Fig. 3. 

The MAE and the orbital magnetic moment were calculated by means of self-consistent 
fully relativistic calculations using the bond lengths obtained in the scalar relativistic cal- 
culations. The MAE was defined as 

MAE = Etot[\\] — Etot[-L] , (2) 

where -Ef t[||] and .E^T] denote total energies of states with magnetization direction parallel 
and perpendicular to the Bz plane, respectively. The choice of the direction parallel to the 



plane is arbitrary, since the in-plane anisotropy is negligible on the scale of the considered 
energies. Results obtained with and without OP corrections are reported. In the former 
case, the spin-dependent OP correction^ was applied to the 3<i-orbitals of Fe, Co, Ni and 
to the 4d-orbitals of Ru, respectively. 

In order to cross-check one of the most important details of the GGA calculations, the 
bonding behavior of the Co dimer with Bz, we also performed ab initio quantum chemical 
calculations at the level of second order M0ller-Plesset perturbation theory (MP2). The 
MP2 results were obtained from the MP2 implementation of Gaussian03.— For C02BZ, the 
Co atoms were described with a scalar-relativistic effective core potential (ECP) replacing 
10 core electrons (MDF10),— with the corresponding (8s7p6d If)/[6s5p3d If] GTO basis 
set of triple-zeta quality. Accordingly, for benzene the Dunning correlation-consistent basis 
sets of double- and triple-zeta quality (cc-pVDZ and cc-pVTZ, respectively)^ were used. 
Since all of these basis sets are rather large, all energies have been corrected for the basis 
set superposition error (BSSE) with respect to the dissociation of the Co dimer from the 
benzene ring, employing the counterpoise scheme proposed by Boys and Bernardt^^ as 
implemented in Gaussian03. 

III. RESULTS AND DISCUSSION 
A. Structure and spin state 

Dimer adsorption energies for each considered total spin and symmetry are shown in 
Fig. [2j The optimized structure parameters, dimer adsorption energies, and related spin 
magnetic moments for the six ground-state structures are listed and compared with literature 
data in Table [D 

1. 3d transition-metal complexes 

We find that the adsorption mode with the dimer axis perpendicular to the benzene 
plane results in the most stable structure for all investigated 3d systems, Fe2Bz, C02BZ, 
and Ni2Bz. This ground-state geometry is consistent with other GGA results obtained for 
dimers on graphene.— The only discrepancy occurs for the ground-state spin magnetic 
moment of Fe2Bz. Here, we found a reduction of the free dimer spin (S = 3) to S = 2, 
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whereas Johll et al. and Cao et al. had reported 5 = 3 for the marginally different situation 
of Fe2 on graphene.— ^ As the related energy difference in our calculation was very small 
(9 meV), we repeated the calculations with full optimization of the C-C and C-H distances 
(Cq v symmetry). The full optimization inverted the order of the two considered states, the 
state with S — 3 now being 16 meV lower than the competing state with 5 = 2. Moreover, 
a ferrimagnetic state with S — 1 is found only about 40 meV higher in energy, Fig. [21 Such 
small energy differences cannot guarantee the stability of the Fe2Bz magnetic ground state 
and should give rise to strong spin fluctuations. 

Turning our attention to the ground-state geometry of Fe2Bz, we note that Parker— pro- 
vided evidence of the same proposed geometry by obtaining an excellent agreement between 
calculated and experimental infrared spectra. Note, that in this geometry Fe2Bz complexes 
with 5 = 1,2, and 3 give almost identical simulated spectra.— Thus, the comparison cannot 
be used to distinguish the magnetic state. Our results add additional weight to Parker's 
arguments, who proposed a Cq v geometry for Fe2Bz. One should note, however, that other 
structure types compete with the Cq v geometry, see Fig. [2j Indeed, if the C and H coordi- 
nates are optimized as well, the ground state turns to ||b, 5 = 3, almost degenerate with 
the states _L C , S = 2 and S = 3. In line with this finding, Cao et al. reported a ground 
state with the Fe dimer above a graphene hollow site, but not perpendicular to the plane.— 
In the following discussion, we will disregard Fe2Bz structures different from Cq v due to the 
mentioned experimental evidence of this structure type. 

For C02BZ and M2BZ, we obtained both the magnetic moment and the ground-state 
structure in agreement with the results by Johll et a/.— and by Cao et al^- The GGA 
dimer adsorption energy, E a d = 1.39 eV for C02 on Bz from our calculation, and the related 
energies 0.92 eV from Johll et al. and 1.13 eV from Cao et al. for C02 on graphene indicate 
a reasonable stability of this structure. Noteworthy, any other spin state considered in 
our calculations, including a ferrimagnetic solution with a total spin 5 = 1, has a much 
higher energy, at least 0.80 eV above the ground state. The ground state of Ni2Bz is also 
sufficiently separated from other states with a different spin and/or geometry, see Fig. [21 
We also checked a further geometry with Ni atoms close to next nearest C-bridges (see Fig. 
Al (iv)). Such a geometry was found by Rao and Jena^ to be lowest in energy. We could 
not confirm this finding and obtained, even with full relaxation, a 0.45 eV higher energy for 
this geometry (5 = 1) than for the _L C geometry with 5 = 1. 
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To investigate the stability of the proposed perpendicular adsorption mode for the Co 
system in more detail, three further structures were considered, including two dissociated 
cases: (i) attachment of one Co atom on each side of the carbon ring, (ii) dissociation of one 
of the Co atoms resulting in one free Co atom and CoBz, and (iii) dissociation of both Co 
atoms resulting in two free Co atoms and one free benzene molecule. In the related scalar 
relativistic atom calculations, non-integer occupation of the open shells was admitted. All 
the three structures have higher energies, by 2.47 eV, 3.49 eV, and 5.05 eV, respectively, 
than the perpendicular arrangement. From the energy difference between state (ii) and state 
(iii), a value of 1.56 eV is found for the adsorption energy of a single Co adatom on benzene. 
We also can find a binding energy of 3.66 eV for the Co dimer from the third dissociated 
state. In comparison with the DFT data, the experimental binding or adsorption energies 
are considerably smaller: 0.34 eV— (adsorption energy of a Co adatom to benzene) and 
1.72 eV— (binding energy of a Co dimer). This is in line with the known tendency of DFT 
calculations to overestimate the binding energies in many cases. 

Thus, to confirm the qualitative validity of the energies and structure sequence, we per- 
formed quantum chemical (MP2) calculations. In comparison with our previously published 
results,— the MP2 calculations were improved by taking into account BSSE corrections and 
by extending the benzene basis. Fig. [3] compares the energy sequences of three possible 
high-symmetry structures and one dissociated configuration obtained by GGA with related 
MP2 results. All MP2 energies were evaluated by single point calculations using the GGA- 
derived geometries, except for the Co dimer. For the latter the interatomic distance was 
optimized at the MP2 level. It turned out slightly shorter (0.1909 nm) than the GGA result 
(0.1997 nm). For all structures, MP2 calculations were carried out for S = 0, 1, 2 and 3. 
In all cases a total spin of 2 was found to be most favorable. Importantly, the adsorption 
energy of the Co dimer to benzene was found to be yet higher than in the GGA calculation. 
It is obvious from Fig. [3] that the quantum chemical calculations confirm the main GGA 
result, that bonding of a Co dimer with a single molecule of benzene results in the structure 
depicted at the bottom of Fig. [3] with a total spin S = 2. Also, the sequence of the higher- 
energy structures is the same in GGA and in MP2, and the same spin magnetic moments 
are found with the exception of the dissociated state, where MP2 predicts S = 3 and GGA 
yields 5 = 2. When improving the benzene basis from a double- to a triple- zeta level, the 
BSSE-corrected adsorption energies rise consistently by about half an eV, reaching 2.36 eV 
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for the most stable structure. This result should serve as a valid proof that the Co dimer 
can be bound to the benzene ring. 

As a final check that the calculations described above provided the correct ground-state 
geometry, we performed a cross-check for C02BZ with the pseudopotential code ESPRESSO- 
4.01.— We carried out a full optimization of all atomic positions starting from 14 kinds of 
initial structures. As before, GGA and a scalar relativistic mode were used. Other technical 
details are described in the appendix. The results confirm that the bonding of C02 with a 
single molecule of benzene very likely results in the perpendicular configuration, which is 
separated from other possible arrangements by at least several hundred meV. Only a tiny 
distortion of the benzene plane is found in the full optimization results. This is a weak Jahn- 
Teller effect that splits the singly occupied two-fold degenerate HOMO state originating from 
C02. This splitting is very small (<1 meV), since the original symmetry of the Co dimer 
where the HOMO resides is only weakly distorted by the hexagonal ligand. Thus, it will 
hardly have any influence on the magnetic properties of C02BZ. In particular, if spin-orbit 
coupling is taken into account, the HOMO will be split by this interaction rather than by 
the Jahn- Teller effect, which is almost 2 orders of magnitude weaker than the spin-orbit 
coupling in the considered case. Only if the magnetization is oriented perpendicular to the 
dimer axis, the spin-orbit splitting vanishes in lowest order. In this case, the Jahn- Teller 
effect might marginally reduce the total energy and, thus, the MAE. 

We conclude that C02BZ and N^Bz probably exhibit a Cq v symmetry, like Fe2Bz. As 
distinct from the Fe system, C02BZ and M2BZ have a stable magnetic ground state. Experi- 
mental evidence of these proposed structures seems however lacking at the moment. For the 
anion Co2Bz~, the observed photoelectron spectra^ allowed to exclude a structure where 
the two Co atoms are placed on both sides of the benzene plane. 

2. 4d transition-metal complexes 

We find (Fig. [2]) that among the investigated Ad dimers only RU2 prefers an upright 
adsorption mode. It binds to the benzene as strongly as the cobalt dimer, but its magnetic 
ground state has a lower spin, S — 1, and lies only 0.25 eV below a zero spin state. The 
ground states of both Rh 2 Bz and Pd2Bz are found to be almost degenerate with respect to 
spin multiplicity (Rl^Bz) or geometry (Pd2Bz). 
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We are not aware of any published information about the geometry of the neutral com- 
plexes RU2BZ and RI12BZ. For the cation Rli2Bz + , the structure type ||b (Fig. [I]) was 
obtained as the lowest-energy structure by DFT calculations using the B3LYP functional.— 
This is the same ground-state structure as we find for the neutral Rh 2 Bz. 

Liittgens et al— deduced the vibration energies of both Pd2Bz and Pd2Bz~ from photo- 
electron detachment spectra. They postulated an orientation of the Pd dimer perpendicular 
to the benzene ring because the observed vibration frequency of Pd2Bz is close to that of the 
free Pd 2 . One should note that our calculated ground-state geometry of Pd 2 Bz contradicts 
this analysis. On the other hand, calculations by Cabria et al— using LSDA and GGA and 
by Loboda et al— using the B3LYP functional find a ground-state geometry with the two 
palladium atoms placed horizontally above the carbon ring, similar to our results. 

We performed a series of additional tests in order to clarify this discrepancy between 
experiment and theory. First, we checked the influence of spin-orbit interaction and found, 
that related total energy shifts do not exceed 0.1 eV. The parallel adsorption mode hence 
is still more stable than the perpendicular one. Second, we checked a possible asymmetric 
adsorption above a single bridge site. Indeed, the related total energy is about 0.5 eV lower 
than for the adsorption above the hollow site, if every symmetry constraint is released. The 
dimer axis then deviates from the initial perpendicular orientation and forms an angle of 
about 40 degrees with the benzene plane, while 5 = 0. Yet, the energy of this structure is 
still higher than that of the parallel configuration. Third, we calculated vibration frequencies 
of the Pd-Pd bond using the harmonic approximation for the three adsorption modes defined 
in Fig. [TJ If 5 = 1, the three structures give rise to almost the same vibration energies of 21.9 
meV (J- c ), 22.3 meV (||b) and 22.5 meV (|| t ). For the asymmetric bridge site configuration, 
a value of 20.6 meV is obtained. All of these energies are close to the vibration energy of 
the free dimer, 26 meV.— This comparison shows that proximity of the vibration spectrum 
to that of the free dimer is no proof of the perpendicular geometry. For the ground state, 5 
= 0, the parallel structures turn out to be much softer, 12.8 meV (||b) and 13.3 meV (|| t ), 
than the free dimer and also than the perpendicular geometry, 19.7 meV. On the whole, 
the calculated vibration spectra do not provide enough evidence in favor of any investigated 
structure. Finally, we optimized the structures of Pd2Bz~ anions. The structure ||b with 5 
= 1/2 is again more stable than the structure _L C , by 0.45 eV. The structure with asymmetric 
bonding above a single bridge site is 0.13 eV lower in energy than the structure _L C , but it is 
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still 0.32 eV higher than the structure ||b- Summarizing this point, the discrepancy between 
experimental and theoretical results on Pd2Bz persists. 

B. Electronic structure and bonding mechanism 

Free TM dimers have been discussed in detail recently.-^^ 1 ^ Analysis of their electronic 
structure reveals that a singly occupied HOMO which is two-fold degenerate in the absence 
of spin-orbit coupling is responsible for the giant magnetic anisotropy predicted in some of 
these dimers.- For example, the most important feature in C02 is a two-fold degenerate singly 
occupied 3d-5* state. It is split by spin-orbit interaction, if the magnetic moment is oriented 
along the dimer axis but stays degenerate if the moment is oriented perpendicular to the 
axis.- Concerning the bonding between metal atoms and a benzene molecule, Mokrousov et 
a/.- reported a schematic analysis for V-Bz complexes and showed that the HOMO and the 
lowest unoccupied molecular orbital (LUMO) of benzene interact with the metal s and d 
orbitals of the same symmetry. 

To better understand the bonding mechanism between the TM dimer and the benzene 
molecule, we compare the levels of C02BZ and of Fe2Bz with those of the related free dimers, 
Fig.|H The third panel (from left) shows the textbook electronic structure of benzene and the 
leftmost panel refers to the free C02, as recently discussed in Ref. fj. The other panels show 
the electronic structure evaluated for the ground-state geometries and spin multiplicities of 
C02BZ, Fe2Bz, and Fe2 as well as for the low-lying S = 3 state of Fe2Bz. It turns out that 
bonding of C02 on benzene does not lead to any deterioration of the magnetic properties of 
C02: (i) the ground-state spin stays S = 2, as in the free dimer; (ii) the Co 3g?-<5* level is 
still two-fold degenerate in C02BZ due to the Cq v symmetry; (iii) this level is still the singly 
occupied HOMO. In this way, the key feature responsible for the giant MAE of the free Co 
dimer is preserved in the C02BZ structure if the benzene molecule binds perpendicularly to 
C02 in Cq v symmetry. Yet, there is an important difference between free C02 and C02BZ. In 
the free dimer, the two Co atoms contribute equal weights to the minority spin 3<i-5* state. 
At variance, the HOMO of C02BZ receives 94% of its weight from that Co atom which is 
farther away from the benzene plane, see Fig. [51 

One can note that the magnetic moment of free C02 is also preserved in the parallel 
arrangements ||b and || t (Fig. [2]). However, the reduction of the symmetry to Civ causes a 
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split of all 5 and ir states (not shown here), resulting in non-degenerate HOMO and LUMO 
in these structures. The minority spin it* state near to the Fermi level is split by 0.57 eV in 
the ||b structure and by 0.60 eV in the || t structure, while the 8* state in the minority spin 
channel is split yet more strongly, by 0.68 eV in ||b and by 0.85 eV in || t . 

We find that, unlike in the case of C02, the adsorption of Fe2 on benzene results in a 
change of both magnetic moment and electron configuration, as compared with the free Fe 
dimer. The ground-state level scheme of Fe2 is very similar to that of C02, see Fig. HI but 
two more holes are introduced in the minority spin channel. As a result, S = 3 (Fe2) instead 
of S = 2 (Co 2 ) and the exchange splitting is enhanced, so that the majority spin Asa* level is 
again quasi-degenerate with the HOMO. The latter is now allocated to the singly occupied 
S orbital instead of 5* in C02. If Fe2 binds to benzene, the electron occupying the majority 
spin it* level of Fe2 moves into the minority spin S level. The minority spin S level becomes 
doubly occupied, while the majority spin n* level turns singly occupied and acts as the 
HOMO in the S = 2 ground state of Fe 2 Bz. 

The electronic structure of Fe2Bz is very similar to that of C02BZ in the same spin state 
(S = 2). The covalent splittings of the d-states are somewhat larger in the Fe system than 
in the Co one, due to the larger extension of the Fe-d orbitals compared with the ti-orbitals 
of Co. This yields a somewhat different orbital order. It will be recalled that in Fe2Bz a 
state with S = 3 is close in energy to the ground state; in this state the exchange splitting is 
larger and both holes enter the minority spin channel, like in the free Fe dimer. We included 
the corresponding level scheme for comparison in Fig. HI The essential difference between 
Fe2 and excited Fe2Bz in the same spin state is that a and S in the minority spin channel 
interchange their positions. Thus, Fe2Bz in the state S = 3 has a fully occupied 5 HOMO 
and, thus, a small magnetic anisotropy. 

The orbital characteristics at the Fermi level determine the main physical properties of 
the structure, while the stability of the complex depends on how the TM dimer binds to the 
benzene molecule. Fig. [5] shows the integrated density of states (IDOS) for the ground state 
of the C02BZ structure and the orbital composition of each state. 

In an axial symmetry, the five d orbitals are split into three groups: da(d z 2), two-fold 
degenerate d7i(d xz , d yz ), and two-fold degenerate dd(d xy , d x 2_ y 2). The six benzene tt orbitals 
can also be classified with respect to the same axis^ and include one La(ni) orbital, two 
degenerate HOMO Lrc{jX2, 7r 3 ), two degenerate LUMO L5(nl, 7r|) , and one L(p{rcl) orbital, 
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where L means ligand. The adsorption between the Co dimer and the benzene molecule is 
realized primarily by forming three types of chemical bonds, 5, it, and sa, between carbon 
atoms and Col (the Co atom nearest to the benzene ring), while the the other Co atom, 
Co2, mainly binds with Col and scarcely contributes to the bonding with the benzene. For 
example, the minority spin 5 state at -4.3 eV and the related 7r*-dominated (Bz-L<5) 5-state 
at -1.2 eV in Fig. [5] stem from the combination of the d xy , d x 2_ y 2 orbitals of Col and the 
LUMO of benzene. The tt states positioned near -7.8 eV in the minority spin channel and 
around -8.0 eV in the majority channel can be attributed to the d xz , d yz orbitals of Col and 
the HOMO of benzene. The hybridization between the 4s orbital of Col and the La(it\) 
orbital of benzene forms the sa orbital, located at around -10.2 eV. Formation of these three 
kinds of bonds lowers the energy of the Co 2 Bz complex compared with the dissociated state 
and stabilizes the perpendicular adsorption structure. 

A simple model for C02BZ can be sketched from above analysis of the bond mechanism: 
the benzene molecule plays the role of a substrate to fix the Co dimer, the Co atom next to 
the benzene plane acts as "glue" to bind the dimer to the benzene, and the other Co atom 
mainly contributes to the states near the Fermi level and dominates the magnetic properties 
of the whole structure. 

Level schemes for the ground states of N^Bz, RU2BZ and for the related free 
dimers/benzene are shown in Fig. [61 The electronic structure of M2 is almost unchanged 
upon adsorption. The spin remains S = 1 and ir* of the minority spin channel is still the 
singly occupied HOMO. The situation is different for RU2, where the total spin is reduced 
from S — 2 to S — 1 due to the adsorption. Together with the change of spin, an essential 
alteration of the orbital order close to the Fermi level takes place. In particular, majority 
spin a* and tt* levels interchange their positions. Thus, RU2BZ has a singly occupied 5* 
HOMO like Co 2 Bz. 

C. Spin and orbital moments 

Site-resolved spin moments fi$ an d orbital moments /ix, for the stable magnetic structures 
with perpendicular geometry (TM = Fe, Co, Ni, Ru) are listed in Table HTl The calculations 
were carried out within the fully relativistic scheme. The moments are aligned perpendicu- 
larly (_L) or parallelly (||) to the benzene plane, respectively. We checked that the in-plane 
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anisotropy of the magnetic moments is marginally small and can be safely neglected. The 
effect of the OP correction is also given for comparison. It is clear that the magnetic moment 
is mainly distributed between the two transition-metal atoms. The magnetic moment on 
the C sites is so small that it can be neglected (|/is(c)| < 0.05/zb). 

Inspection of the fis data reveals that in all four systems the TM2 atom shows a much 
higher spin moment than the TM1 atom. The relatively high coordination number of TM1, 
seven, results in a considerable reduction of its spin moment. On the other hand, TM2 is 
only singly coordinated and thus behaves almost like a free atom. The spin moments of Fe, 
Co, Ni, and Ru free atoms amount to 4, 3, 2, and A fig, respectively. The calculated spin 
moments of the TM2 atoms in 3<i-TM 2 Bz complexes are only 0.5. . .0.7 fiB smaller than 
the corresponding atomic values. In the case of Ru, the TM2 carries only about half of 
the atomic spin moment. Spin magnetism of Ad atoms is in general less stable than that of 
isoelectronic 3d atoms, since the Ad intra-atomic exchange (Stoner) integrals are somewhat 
smaller than the related 3c? integrals. 

Another feature is that the spin moments are nearly the same for both magnetization 
directions. When the moment orientation switches from _L to ||, the primary effect is a 
change of the orbital moment of the TM atoms. This fact indicates that the magnetic 
anisotropy of these systems is closely connected to the anisotropy of the orbital moments.— 
It is worth noting that in both the C02BZ and the RU2BZ systems the TM2 atoms show very 
large orbital moments in the _L orientation and relatively small values in the || orientation. 
This is a sign of a large magnetic anisotropy of these systems. 

As expected, the calculated moments in Table [III show that the OP correction generally 
increases the orbital moment while scarcely affecting the spin moment. For example, when 
the OP correction is allowed for, the orbital moments of Co atoms with magnetic moments 
parallel to the Bz plane are about three times larger than those calculated without the OP 
correction. This is caused by the very construction of the OP correction scheme, where 
additional (exchange) energy is gained if the orbital moment is enhanced.— 1 ^ In case of 
perpendicular orientation of the moments with respect to the Bz plane, the orbital moments 
are less influenced by the OP correction, since spin-orbit coupling alone already provides 
almost the maximum orbital moment allowed by the given electronic level sequence. 

The total orbital moments evaluated for the case when the magnetic moment is parallel 
to the dimer axis (i.e., perpendicular to the benzene plane) directly reflect the nature of the 
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HOMO. In the case of Fe2Bz, Fig. HJ the HOMO is a tt* state in the majority spin channel. 
Spin-orbit coupling splits this state in such a way that the energy of the m — — 1 sub- 
level is reduced (m denotes the magnetic quantum number). This sub- level is consequently 
occupied, while the m — +1 sub- level is empty, and the total orbital moment is close to 
— 1 fiB- in the case of C02BZ, Fig. HJ the HOMO is a 5* state in the minority spin channel. 
Spin-orbit coupling splits this state so that the energy of the m = +2 sub-level is reduced. 
Accordingly, [il ~ 2 \Ib- In M2BZ, the HOMO is a 5* state with both sub- levels occupied. 
Therefore, the orbital moment nearly vanishes. Finally, in RU2BZ the HOMO is a 5 state in 
the minority spin channel with fi L £s 2 ji B . 



D. MAE 



After analyzing the spin and the orbital moments in the stable _L C TM2BZ structures, we 
proceed to another important property of magnetic systems, the MAE. This quantity is in 
the main focus of the present investigation. In the perpendicular adsorption mode of TM 2 
on the benzene molecule, it is natural to consider the MAE as the energy difference between 
the states with magnetization direction parallel and perpendicular to the benzene plane, Eq. 
©. 

The first two lines in Table II III list our calculated MAE for the stable magnetic TM2BZ 
structures with TM = Fe, Co, Ni, Ru. The data for E to t were obtained by two self-consistent 
fully relativistic calculations with respective magnetization directions. The values calculated 
without OP correction should be considered as a lower estimate to the expected MAE, while 
an upper estimate is obtained by including the OP correction. 

Endowed with a large ground-state orbital moment as demonstrated in Sec. IIII C[ C02BZ 
and RU2BZ show a huge MAE. The lower estimate to the MAE in C02BZ is hardly changed in 
comparison with the free Co dimer.— This is because the magnetic state and the important 
features of the electronic structure of C02 are not changed by the adsorption. The upper 
estimate, 334 meV per C02BZ molecule, is even higher than the related value for Co 2 (188 
meV per dimer). This is due to the almost complete localization of the HOMO on TM2 in the 
case of C02BZ (Table HTT1 lines three and four). While /^(tm2) ~ 2 /xb in C02BZ (Table HT1). it 
is only half as large in the free Co dimer, where by symmetry /iL(TMi) = Hl(tm2) ~ 1 /is in 
the ground state. The OP correction energy is only half as large in the latter case compared 
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to the former one, since it is quadratic in the atom-projected //£,. At variance, the spin-orbit 
coupling energy is linear in ^l- 

The case of Ra^Bz is different. We find in the ground state of the free Ru dimer (with S 
= 2) a two-fold degenerate, completely occupied majority-spin 3d — 7r* HOMO and, thus, a 
small MAE. The spin moment is reduced by the adsorption of RU2 on benzene. This causes 
a change of the electron configuration, resulting in a two-fold degenerate and singly occupied 
5* HOMO with a related huge MAE. 

The iron system shows a somewhat smaller MAE, mainly due to the smaller value of \m\ 
of the HOMO. Finally, the Ni system has a fully occupied HOMO that does not contribute 
to the MAE in lowest order. Nonetheless, the obtained MAE reaches or exceeds the highest 
known experimental values.- 1 ^ 

We finally would like to understand the obtained numbers in terms of simple arguments 
based on perturbation theory,- extended here to include the OP correction. Our consider- 
ation is limited to systems with singly occupied, two-fold degenerate HOMO (here, Fe2Bz, 
C02BZ, and Ru2Bz). The MAE is approximated by the single-particle energy change of 
the occupied HOMO level upon changing the direction of magnetization. If the OP correc- 
tion is included, we call the magnetic anisotropy energy MAE^o+op); otherwise it is called 
MAE ( so). Thus, 

MAE ( so) « \m\ £ (C^_ TMi + C 2 _ m4 _ TM j &/2 (3) 

i 

in first-order perturbation theory. Here, ^ is the <i-shell spin-orbit parameter and Cm^-TMi 
is the projection of one of the HOMO orbitals on the <i-orbital of atom TMz (i = 1, 2) with 
magnetic quantum number m (\m\ = 1 or 2 for a HOMO of type it or 5, respectively). If 
OP corrections are taken into account, this first-order estimate changes to 

MAE (so+ op) ~ H ]T {C 2 m4 

-TMi t ^-m,d-TUiJ 

(U/2 + BA\fj, L(TMi) \/fj, B ) . (4) 

i 

Here, B denotes the TM-specific Racah parameters,— evaluated from the related atomic 
orbitals, and A|/i L{TM i)| = |j"f( T Mi)l ~ l/4,(TMi)l> (* = 2 )> according to Tab. OH 

Table mH shows the major contributions to the composition of the HOMO, (C^ d _ TMi + 
C^ m(J _TMi)j the related magnetic quantum number \m\, the occupation of the HOMO, the 
spin-orbit parameters, and the Racah parameters. MAE values estimated by first-order 
perturbation theory are given in brackets following the self-consistently evaluated MAE 
data. We find that the self-consistent values are smaller than the estimates obtained by 
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perturbation theory. This has at least two reasons: (i) negative contributions from levels 
other than the HOMO and higher-order HOMO contributions, cf. the results for Ni2Bz 
where the HOMO does not contribute in first order, and (ii) charge relaxation reduces 
the effect. Nonetheless, the self-consistent MAE amounts to 60 ... 80% of the first-order 
estimates. 

According to the above analysis, C02BZ and RU2BZ are interesting candidates for strong- 
MAE applications. It should be noted that in both cases the easy axis of magnetization 
is directed perpendicularly to the benzene ring. This fact is advantageous for conventional 
recording techniques. 

IV. SUMMARY 

We report a systematic DFT study of the ground-state structures, bonding mechanism, 
spin and orbital moments, and in particular of the MAE of TM2BZ complexes (TM=Fe, Co, 
Ni, Ru, Rh, Pd), using the full-potential local-orbital method FPLO. Upright adsorption 
modes with Cq v symmetry of Fe2, C02, N12, and RU2 on benzene molecules are confirmed 
(Fe2Bz) or predicted (others). Huge MAE, stable geometry, and stable magnetic ground 
states are predicted for C02BZ and for RU2BZ. The main origin of the large anisotropy of 
these two systems is the large orbital moment of the TM atom which is farther away from 
the benzene plane. Analysis of the electronic states shows that bonding of the Co dimer on 
the benzene molecule does not lead to any deterioration of the magnetic properties of the 
dimer. Most important is that the two-fold degenerate singly occupied HOMO state of the 
free dimer is preserved, which allows the spin-orbit coupling to produce a large magnetic 
anisotropy. An important conclusion can be drawn from these results: robust and easy-to- 
prepare carbon-based substrates are well-suited to adsorb transition-metal dimers for the 
purpose of high-density magnetic recording. We hope that our predicted exceptionally large 
MAE of C02BZ and RU2BZ will motivate experimental investigations of transition-metal 
dimers on carbon-based substrates, like graphite or graphene. 
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Appendix 



A full optimization of all atomic positions was carried out for C02BZ complexes by us- 
ing the ESPRESSO code.— We used the pseudopotentials Co.pbe-nd-rrkjus.UPF, C.pbe- 



rrkjus.UPF and H.pbe-van_bm.UPF from the http://www.quantum-espresso.org distri- 
bution. A supercell of the size 20 Ax 20 Ax 20 A was used to make sure that there is 
virtually no interaction between the molecules of neighboring cells. The Brillouin zone sam- 
pling was performed only on the T point. The cutoffs used for the wave functions and for 
the charge density were 60 Ry and 300 Ry, respectively. The convergence in total energy 
was carefully checked. A Marzari-Vanderbilt cold smearing^ with 0.007 Ry was used to 
get the convergence in energy levels. The optimizations were done without any constraints 
on symmetry or spin moment. Fig. IA1I shows the 14 initial structures used to start the 
geometry optimization. 

The optimization results in six types of final structures. Structures (i), (xi), (xii), and 
(xiii) converge to the perpendicular adsorption mode; structure (x) stays almost unchanged 
and leads to the parallel adsorption mode on the bridge site of the carbon ring; structures 
(v) and (ix) converge to the parallel adsorption mode on the top site of the carbon ring; 
structures (iii), (vi), and (viii) finally go over to a structure type similar to (iii); structures 
(ii) and (iv) become the structure type (iv); and the final structure of (xiv) is still the 
adsorption mode on both sides. The optimization for structure (vii) does not converge. 
Among these structures, the perpendicular adsorption mode shows the lowest total energy, 
which is 0.74 eV lower than the ||b adsorption mode. The spin of this structure is S = 2, 
confirming the FPLO result. There, the ||b adsorption mode was found 0.98 eV higher than 
the perpendicular adsorption mode. 

Full optimization of all atomic positions shows that there is only a tiny distortion of the 
benzene plane in the perpendicular adsorption mode. The length of the C-C bonds increases 
slightly, two of them changing from 1.40 A to 1.4175 A and four of them changing to 1.4169 
A. The C-H bond-length changes from 1.09 A to 1.0889 A and 1.0887 A. The electronic 
structure shows that the two 5* states (now, HOMO and LUMO) are still nearly degenerate 
with a gap smaller than 1 meV. 
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(a) _L C (b) || b 

FIG. 1: (Color online) Three possible high-symmetry structures for TM2BZ complexes. Hexagons 
and blue bullets indicate benzene rings and and transition-metal atoms, respectively, (a) _L C — the 
TM dimer is situated on the C% v symmetry axis perpendicularly to the benzene plane; (b) ||b — 
the TM dimer is parallel to the benzene plane with TM atoms near the middle of opposite C-C 
bonds [p2v symmetry); (c) ||t — the TM dimer is parallel to the benzene plane with TM atoms 
near the top of opposite C atoms {p2v symmetry). 



(c) lit 
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fixed total spin moment [\x ] 



FIG. 2: (Color online) Scalar-relativistic dimer adsorption energies E a d calculated for optimized 
structures of (a) Fe2Bz, (b) C02BZ, (c) M2BZ, (d) RU2BZ, (e) Rl^Bz, and (f) Pd2Bz complexes. 
For all systems, three initial structures illustrated in Fig. [1] were optimized with fixed C and H 
coordinates for the following values of the fixed total spin moment, ns=0, 2, 4, and 6^b- Both 
parallel and anti-parallel relative spin orientations were considered. Open (filled) symbols indicate 
that the moments of the two transition-metal atoms in the spin state with the highest adsorption 
energy are parallel (anti-parallel). 
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FIG. 3: (Color online) Energies and magnetic moments of different C02BZ configurations calculated 
by DFT (left column) and by MP2 (right columns). DZ and TZ abbreviate double-zeta and triple- 
zeta basis sets, respectively. The energies refer to the ground-state energy. For the uppermost, 
dissociated and for the lowermost, ground-state configuration, the spin moments in red (on the 
right-hand side of the Co dimer) refer to the MP2 calculations and the moments in black (on the 
left-hand side of the dimer) refer to the GGA calculations. For the other configurations, GGA and 
MP2 yield the same spin. Hexagons and blue bullets indicate Bz and Co, respectively. 
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FIG. 4: (Color online) Scalar-relativistic single-particle levels of C02 (left panel), C02BZ (ground- 
state structure, second panel), benzene (third panel), Fe2Bz (ground-state structure, fourth panel 
and first spin-excited state, fifth panel), and Fe2 (right panel). All energies refer to a common 
vacuum level. Black lines denote occupied states, orange (gray) lines denote empty states, and 
thick blue (gray) lines indicate singly occupied two-fold degenerate states. With the exception 
of benzene, the levels are spin-split (S = 2, 2, 2(3), and 3 for C02, C02BZ, Fe2Bz, and Fe2, 
respectively). Majority states are indicated by up-arrows, minority states by down-arrows. Dimer- 
dominated states are labeled in black and benzene-dominated states are labeled in red (gray). 
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FIG. 5: The orbital composition of each state (left axis) and the integrated density of states 
(IDOS, right axis) for the ground-state structure and spin of C02BZ. The upward (downward) 
arrow indicates majority (minority) spin states. All energies refer to a common vacuum level. The 
labels for each state are the same as those in Fig. UJ The three types of chemical bonds between the 
Co dimer and benzene are labelled blue in parentheses. The Co atom closer to the benzene plane 
is labelled Col, the other one is labelled Co2. The position of the Fermi level (Ep) is indicated by 
a vertical line. Missing contributions to the orbital composition which should add up to unity are 
due to the omitted C-2s and H-states. 
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FIG. 6: (Color online) Scalar-relativistic single-particle levels of M2 (left panel), N12BZ (ground- 
state structure, second panel), benzene (third panel), RU2BZ (ground-state structure, fourth panel), 
and RU2 (right panel). All energies refer to a common vacuum level. Black lines denote occupied 
states, orange (gray) lines denote empty states, and thick blue (gray) lines indicate singly occupied 
two-fold degenerate states. With the exception of benzene, the levels are spin-split (S = 1 in all 
cases but Ru2, where S = 2). Majority states are indicated by up-arrows, minority states by down- 
arrows. Dimer-dominated states are labeled in black and benzene-dominated states are labeled in 
red (gray). 
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FIG. Al: (Color online) Illustration of 14 initial structures optimized by the ESPRESSO code. 
From (i) to (x), top- views of possible parallel adsorption modes; from (xi) to (xiii), side-views of 
possible upright adsorption modes are shown; (xiv) is a case in which one Co atom is attached on 
each side of the carbon ring. 
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TABLE I: Dimer adsorption energies, E a( i, total and atom-resolved spin magnetic moments, 
A*S(totai)> MS'(tmi)) an d / i s , (TM2) (TM1 refers to the atom closer to the benzene in case of per- 
pendicular bonding, TM2 to the other atom), the distance between the two transition-metal atoms 
cZtm-tMi and the distance between benzene plane and TM1, c^tm-Bz, for the ground-state struc- 
tures of TM 2 Bz (TM = Fe, Co, Ni, Ru, Rh, Pd) complexes. C-C and C-H bond lengths are fixed 
(1.40 A and 1.09 A, respectively). The structure type of each molecule is labelled according to the 
notation introduced in Fig. [TJ Our present results are labelled "a" , literature data (for dimers on 



graphene) are labelled "b", Ref. 



40 and "c", Ref. 
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system 
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lb 


lb 


lb 


E ad (eV) 


0.87 


0.72 


1.39 


0.92 


1.12 


0.96 


1.40 


0.97 


1.16 


1.28 


/^S(total) (Ms) 


4 


6 


4 


4 


2 


2 


2 


2 








Ms(tmi)(mb) 


0.75 


2.76 


1.64 


1.66 


0.71 


0.73 


-0.10 


1.01 


0.00 


0.00 


^S(TM2)(^b) 


3.35 


3.48 


2.45 


2.43 


1.29 


1.29 


2.13 


1.01 


0.00 


0.00 


^tm-tm(A) 


2.04 


2.08 


2.09 


2.03 


2.15 


2.14 


2.22 


2.49 


2.80 


2.75 


dTM-Bz(A) 


1.60 


1.86 


1.66 


1.72 


1.71 


1.73 


1.77 


2.08 


2.14 


2.15 
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TABLE II: Spin moments fis and orbital moments hl (in for the ground-state structures of 
TM2BZ (TM = Fe, Co, Ni, Ru) complexes calculated within the fully relativistic scheme with mag- 
netization perpendicular (_L) or parallel (||) to the benzene plane. The effect of the OP correction 
is also illustrated by comparing the values calculated without the OP correction (SO) and with the 
OP correction (SO+OP). 



Fe 2 Bz Co 2 Bz Ni 2 Bz Ru 2 Bz 





SO 


SO+OP 


SO 


SO+OP 


SO 


SO+OP 


SO 


SO+OP 


II 1 - 

^S(TMl) 


0.75 


0.73 


1.64 


1.64 


0.71 


0.71 


-0.09 


-0.10 


^S(TMl) 


0.75 


0.75 


1.64 


1.64 


0.71 


0.67 


-0.09 


-0.09 


^L(TMl) 


-0.72 


-0.89 


0.07 


-0.10 


0.01 


0.02 


0.04 


-0.04 


^l(TMl) 


0.02 


0.05 


0.13 


0.39 


0.04 


-0.31 


0.03 


0.02 


A i 5(TM2) 


3.34 


3.36 


2.45 


2.46 


1.28 


1.28 


2.13 


2.13 


/ i 5(TM2) 


3.34 


3.34 


2.45 


2.45 


1.28 


1.33 


2.10 


2.11 


A i L(TM2) 


-0.17 


-0.17 


1.93 


2.12 


0.01 


0.00 


1.91 


2.01 


/ i l(TM2) 


0.11 


0.25 


0.17 


0.54 


0.37 


1.69 


0.11 


0.21 


/^(total) 


-0.89 


-1.06 


2.00 


2.02 


0.02 


0.02 


1.95 


1.97 


(total) 


0.13 


0.30 


0.30 


0.93 


0.41 


1.38 


0.14 


0.23 
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TABLE III: The MAE (per molecule), calculated using Eq. ([2]) for the ground-state structures of 
TM2BZ (TM = Fe, Co, Ni, Ru). Positive values of MAE indicate that the easy axis of the system 
is perpendicular to the benzene plane, while negative values mean that the direction parallel to the 
benzene plane is the easy axis. Both a lower estimate of the MAE calculated without OP correction 
(MAE(so)) an d an upper estimate of the MAE obtained with OP correction (MAE(go+op)) are 
listed. Data in brackets indicate estimates obtained by first-order perturbation theory, see text. 
Further, the principal composition ^_TMi 

+ Ci 

md-TMi' the magnetic quantum number |m|, 
and the occupation of the HOMO are given as well as the spin-orbit coupling parameter £^ and 
the Racah parameter B. 





Fe 2 Bz 


Co 2 Bz 


Ni 2 Bz 


Ru 2 Bz 


MAE ( so) (meV) 


+15 [+25] 


+51 [+74] 


-7 


+104 [+123] 


MAE (so+ op) (meV) 


+61 [+107] 


+334 [+519] 


-96 


+279 [+403] 


r 2 + r 2 


68% 


3% 


21% 


2% 


r 2 a- r 2 


15% 


94% 


77% 


94% 


|m| of the HOMO 


1 


2 


2 


2 


occupation of the HOMO 


1 


1 


2 


1 


id (meV) 


61 


76 


96 


128 


B (meV) 


140 


149 


154 


95 
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